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We perform a first principles band calculation for quasi-two-dimensional organic superconductors 
^-(BDA-TTP)2l3 and /3-(BDA-TTP)2SbF6. The first principles band structures between the I3 and 
SbFe salts are apparently different. We construct a tight-binding model for each material which 
accurately reproduces the first principles band structure. The obtained transfer energies give the 
differences such as (i) larger dimerization in the I3 salt than the SbFe salt, and (ii) different signs and 
directions of the inter-stacking transfer energies. To decompose the origin of the difference into the 
dimerization and the inter-stacking transfer energies, we adopt a simplified model by eliminating the 
dimerization effect and extract the difference caused by the inter-stacking transfer energies. From 
the analysis using the simplified model, we find that the difference of the band structure comes 
mainly from the strength of dimerization. To compare the strength of the electron correlation 
having roots in the band structure, we calculate the physical properties originated from the effect of 
the electron correlation such as the spin susceptibility applying two particle self-consistent (TPSC) 
method. We find that the maximum value of the spin susceptibility of the I3 salt is larger than that 
of the SbFe salt. Hypothetically decreasing the dimerization within the model of the I3 salt, the 
spin susceptibility takes almost the same value as that of the SbFe salt for the same magnitude of 
the dimerization. We expect that the different ground state between the I3 and SbFe salt mainly 
comes from the strength of the dimerization which is apparently masked in the band calculation 
along a particular /c-path. 


I. INTRODUCTION 

There have been attempts to synthesize strongly cor¬ 
related electron systems in organic conductors by ap¬ 
plying chemical modification to stable metallic donor 
molecules. For example, there are (S',5')-DMBEDT- 
TTF— and meso-DMBEDT-TTF—, where two methyl 
groups are attached to BEDT-TTE, and they are 
pressure-induced superconductors. In the present ar¬ 
ticle, we theoretically study superconductors based on 
BDA-TTP molecule, which is extended to six-membered- 
ring from five-membered-ring in the a-bond framework 
of BDH-TTP molecule^. The actual materials are 13- 
(BDA-TTP) 2 l 3 and (BDA-TTP) 2 SbEe, which will be 
abbreviated as I3 and SbEe salts, respectively. In both 
materials, conductive layer is the BDA-TTP layer, and 
the anion layer separates the adjacent conductive layers 
as shown in Eig. [1] (a). Molecular configuration in the 
conductive layer is the yd-type as shown in Eig. [1] (b). 
Both materials consist of the stacking structure of the 
BDA-TTP molecules. However, they are somewhat dif¬ 
ferent in that the inter-stacking direction is slightly tilted 
in the I3 salt, but almost side-by-side for the SbEe salt, 
which will be shown later. 

The I3 salt is an insulator at ambient pressure, and 
the superconductivity appears around 10 K under hy¬ 
drostatic pressure of above 10 kbar— . Recently, applying 



FIG. 1 : (color online) Crystal structure of the I3 salt from (a) 
the side view and (b) the conductive layer of the BDA-TTP 
molecules. 

uniaxial strain along the c-axis has given higher Ap¬ 
plying the uniaxial compression once increases the Tc and 
takes a maximum before it decreases^. It is considered 
that applying the pressure in the I3 salt increases the 
overlap between the upper and lower bands, which grad¬ 
ually changes the character of the system from a strongly 
correlated half-filled system to a moderately correlated 
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quarter filled system. The c-axis strain more efficiently 
increases the band-width of the overlap. As the electron 
correlation is reduced to some extent by pressure, the 
insulating nature of the material is lost, and supercon¬ 
ductivity appears^. Theoretically, Nonoyama et al. have 
studied the nature of the charge ordering state and the 
pairing mechanisms in the model of the I 3 salt derived 
from the extended Hiickel band structure^. 

The SbFe salt exhibits superconductivity at 7.5 K at 
ambient pressure^. As for the SbFe salt, there have 
been some controversies regarding both the anisotropy 
of the Fermi surface^ii^ and the directions of the nodes 
in the superconducting gapii“— . In our previous studyi^ 
for the / 3 -(BDA-TTP) 2 MF 6 (M=P, As, Sb and Ta), we 
have obtained the band structure from the first principles 
band calculation, and suggested the origin of the differ¬ 
ences from the extended Hiickel band structurei^. Also, 
there have been some studies on pairing mechanisms me¬ 
diated by spin and/or charge fluctuations in the model of 
^-(BDA-TTP) 2 A:. As for the MFq (M= As, Sb) salts, 
adopting models derived from the extended Hiickel calcu¬ 
lation, Nonoyama et al^ have applied random phase ap¬ 
proximation (RPA) to the two band model, while Suzuki 
et al^ have applied the fluctuation exchange (FLEX) 
approximation to the original two-band model and the 
single-band dimer model. Recently, we have constructed 
the tight-binding model derived from the first princi¬ 
ples band calculation, studied the pairing symmetry of 
the gap function within the spin fluctuation mediated 
pairingi^. 

In the present study, given the difference in the ground 
state between the I 3 salt and the SbFe salt, we focus on 
the difference in the electronic structure between the two 
salts. In fact, despite the similar lattice structure, the 
band structure of the I 3 sah4 and that of the SbFe sah4 
obtained by the extended Hiickel method are known to be 
very different. Here, we perform the first principles band 
calculation for / 3 -(BDA-TTP) 2 l 3 and construct an effec¬ 
tive tight-binding model that reproduces the first prin¬ 
ciples band structure. We compare the band structure 
of the I 3 salt to that of the SbFe salt obtained in our 
previous studyi^, and pin down the origin of the appar¬ 
ently large differences. In particular, we study the re¬ 
lation between the strength of the electron correlation 
and the molecular dimerization. We consider the Hub¬ 
bard model by introducing repulsive interaction between 
the electrons on the same BDA-TTP molecule. Then, 
we study the effect of the electron correlation by apply¬ 
ing the two particle self-consistent (TPSC) method, and 
present quantities such as the spin susceptibility against 
the temperature and dimerization strength, which reflect 
physical properties originating from the electron correla¬ 
tion. We conclude that the ground state of the I 3 salt 
differs from that of the SbFe salt due to the strength of 
the dimerization. 


II. METHOD 

A. first principles band calculation and model 
construction 

We perform first principles band calculation using all¬ 
electron full potential linearized augmented plane-wave 
(LAPW) - 1 - local orbitals (lo) method within the frame¬ 
work of WIEN2k2^. This implements the density func¬ 
tional theory (DET) with different possible approxima¬ 
tion for the exchange correlation potentials. The ex¬ 
change correlation potential is calculated using the gen¬ 
eralized gradient approximation (GGA). 

The single-particle wave functions in the interstitial 
region are expanded by plane waves with a cut-off of 
^MxA’niax = 3.0 duc to the presence of the hydrogen 
atom, where Rmt denotes the smallest muffin-tin radius 
and Amax is the maximum value of K vector in the plane 
wave expansion. In the I 3 salt, the muffin-tin radii are 
assumed to be 2.50, 1.62, 1.15, and 0.62 atomic units 
(a.u.) for I, S, G, and H, respectively. Amax is taken as 
4.8, and the plane wave cutoff energy is 318.6 eV. In the 
SbEe salt, the muffin-tin radii are assumed to be 1.74, 
1.74, 1.62, 0.83, and 0.45 a.u. for Sb, E, S, C, and H, 
respectively. Amax is taken as 6.7, and the plane wave 
cutoff energy is 604.7 eV. Calculations were performed 
using 6x3x9 /c-points for the I 3 salt and 7x3x9 /c-points 
for the SbEe salt in the irreducible Brillouin zone. We 
adopt the lattice structure determined experimentally for 
each materials^i^, and we do not relax the atomic posi¬ 
tions in the calculation. 

Having done the first principles band calculation, we 
then construct a tight-binding model which accurately 
reproduces the first principles band structure. Erom the 
lattice structure of the two materials, we regard one 
molecule as a site and consider a two-band (two sites 
per unit cell) tight-binding model to fit the first princi¬ 
ples band structure. The tight-binding Hamiltonian, Hq, 
is written in the form 

Ho = \j^ia:j^cl^^Cjj3cr +H.C.| , (l) 

{ia:jP),a 

where i and j are unit cell indices, a and [3 specifies the 
sites in a unit cell, (ciaa ) is a creation (annihilation) 
operator with spin a at site a in the i-th unit cell, tia-j^ 
is the electron transfer energy between (i,(a) site and 
(j, yd) site, and {ia : jP) represents the summation over 
the bonds corresponding to the transfer. 

By Eourier transformation, eq. (HD is rewritten as 

Ho = ^ (^) 

k ,(T,a,/3 

where 5^/3 {k) is the site-indexed kinetic energy repre¬ 
sented in /g-space. The band dispersion is given by diag¬ 
onalizing the matrix 5^/3 (/s), 

^a/3 w = E daj (k) ik)^Ak), (3) 
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where {k) gives the band dispersion of the 7 -th band 
measured from the chemical potential, and (k) is the 
unitary matrix that gives the unitary transformation. 

We adopt the two-band Hubbard model obtained by 
adding the on-site (intra-molecule) repulsive interaction 
to the tight-binding model derived from the fitting of the 
first principles band structure. The Hubbard Hamilto¬ 
nian, iJ, is 

H = Ho + '^ Uoniatriiai (4) 

ia 

where Uq is the bare on-site interaction and riiaa is the 
number operator of the electron on the a-site in the i- 
th unit cell. Since both salts are configured as a form 
of D 2 X where D is the donor molecule and X~^ is the 
anion, the band-filling is 1/4-filled in the hole represen¬ 
tation (3/4-filled in the electron representation). 


B. Two particle self consistent method 


To deal with the electron correlation effect arising from 
the on-site repulsion, we apply TPSC to the multi-site 
Hubbard model given by eq. as follows. The bare 
susceptibility in the site-representation is given by 

= -^^Gljs{k + q)G%{k), (5) 

^ k 


where T and Nc are the temperature and the total num¬ 
ber of unit cells, respectively, and (k) is the bare 
Green’s function given as 

G%{k) = ( 6 ) 

^ [KJ 


Here, we introduce the abbreviations k = {k^iSn) and 
q = for the fermionic and bosonic Matsubara 

frequencies. The indices ajd means {a yd)-element of the 
matrix such as {q). 

TPSC has been applied to single-site system s^^i^^ , 
multi-site system^ and multi-orbital systen>^. By ap¬ 
plying TPSC, we can consider the local vertex correction 
in both spin and charge channels within a self-consistent 
procedure. In the TPSC, using the bare susceptibility 
given by eq. m, the spin and charge susceptibilities are 
obtained as 


r^{q) = \°{q), (7) 


, 1 -1 


V7?) = i + x°{q)u'^'^ 7(9), 


(8) 


where is the local spin (charge) vertex and I 

is the unit matrix. The local vertices are determined by 
satisfying two sum rules for the local moment such as 

■jr Y. = na-2 (nafriai ), (9) 

= na + 2 {nafUal) - nl, ( 10 ) 

" 9 


where is the particle number at the site a. We have 
used the relations = n /2 and n^a = '^aa from 

the Pauli principles. 

The local spin vertex is related with the double 
occupancy {ria^nai) by the following ansatz 


= 

^ OLOL 


rjO 


( 11 ) 


where is the {a (a)-element of the on-site interaction 
matrix U^. Equation (pT]) breaks the particle-hole sym¬ 
metry and should be used for < 1. When > 1, that 
can be applied through the particle-hole transformation, 
then the double occupancy given by 


TTSp 2 / TTSp \ 

Da = j^^+(l-'^){na-l)e{na-l), (12) 

^aa V ^aa J 

where d (x) is Heaviside step function. Equations (|7|)- (pT]) 
give a set of the self-consistent equations for the TPSC 
method. Obtaining the l/gp and l/ch, the interaction for 
the self-energy is obtained as 


(q) = 2 (q) u° + ( 9 ) u°\ • (i3) 

Using the eq. m, the self-energy is given by 

Sa /3 (7 = T ^ 1/7 (q) Gap (k-q), (14) 


and the dressed Green’s function is obtained as 
G{k) = G° {k) + G° {k)t{k)G{k). 


(15) 


Since we need two sites per unit cell,^ 

7 ( 9 ), V’'(9), X^^iq), V^{q), S(fc), G°{k) and G {k) 
all become 2x2 matrices. In the present study, the spin 
susceptibility is obtained as the larger eigenvalue of the 
2x2 spin susceptibility matrix. We consider not only the 
spin susceptibility, but also other physical values such as 
the local spin vertex and the double occupancy. In the 
present calculation, we take the system size as 64 x 64 
k-meshes and 16384 Matsubara frequencies. 


III. RESULTS 

A. first principles band calculation 

Eigures [ 2 ] (a) and (c) show the first principles band 
structures for the I3 and the SbEe salts. Eor both mate¬ 
rials, the experimental lattice structure at an ambient 
pressure and room temperature are used. In both of 
the materials, it can be seen that the highest-occupied 
molecular orbital (HOMO) is isolated from the lowest- 
unoccupied molecular orbital (LUMO). Considering this 
and also the number of donor molecules in a unit cell, 
we adopt the HOMO and HOMO —1 bands as the target 
bands to construct an effective tight-binding model. 
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FIG. 2: (color online) (a) Calculated first principles band 
structure and (b) Fermi surface for the I 3 salt, (c) first prin¬ 
ciples band and (d) Fermi surface for the SbFe salt. In both 
figures of the band structures, the red curves represent the 
first principles bands and the blue solid curves gives the tight- 
binding fit. 


Although the difference is only the anion, the band 
structures of the two materials are apparently very differ¬ 
ent. In order to reveal the origin of this difference in the 
band structure, in the following we focus on the following 
two differences of the two salts. One is the magnitude of 
the molecular dimerization, namely the dimerization of 
the donor molecule in the I3 salt is larger than that in the 
SbFe salt resulting in a larger gap between HOMO and 
HOMO —1 in the former. The other is the anisotropy of 
the band structure, namely, there are two flat portions 
near the Fermi level around the Z and the X-points in 
the I3 salt, while there is only one fiat portion around 
the B-point in the SbFe salt. 

Figure [ 2 ] (b) shows the Fermi surface of the first prin¬ 
ciples band calculation for the I3 salt, where the high 
symmetry points in the Brillouin zone are presented only 
on the kyikij) = 0 plane. The Fermi surface of the I3 salt 
is disconnected, namely quasi-one-dimensional, but it is 
actually close to two dimensional because a slight shift 
of the band structure around the Z-point would give a 
closed (i.e. 2 D) Fermi surface. Figure [ 2 ] (d) shows the 
Fermi surface of the SbFe salt. The Fermi surface is 
cylindrical, reflecting the two-dimensionality of this salt 
as shown in our previous worki^. 


B. Effective tight-binding model 

Figure [ 3 ] shows the effective tight-binding model 
adopted to fit the first principles band. The nearest- 
neighbor transfers are shown in the left panel of Fig. 
[ 3 l and in addition we also need to introduce the next- 
nearest-neighbor transfers shown in the right panel of 
Fig. [ 3 ] to reproduce the first principles band structure 
more accurately. Note that the stacking direction of the 
BDA-TTP molecules is taken in the a-direction^i. The 
band dispersions of the tight-binding model are shown as 
blue solid curves in Fig. [ 2 ] (a) for the I3 salt and Fig. [ 2 ] 
(c) for the SbFe salt. 


^ ^ ^ ^ ^ ^ 



FIG. 3: (color online) The tight-binding model for /3-(BDA- 
TTP) 2 X, where left (right) panel shows the first (second) 
nearest-neighbor transfer energies. 

The transfer energies for the two salts are summarized 
in Table m The bottom three lines represent the magni¬ 
tude of the dimerization which is measured by the ratio 
tp2ltpi^ and the transfer between the inter-stacking di¬ 
rection normalized by the average value of intra-stacking 
transfer energies, (tgi+tg2)/(^pi+^^2) and 2tc/(tpi+^^2)- 
From Table HI it can be seen that there are two major 
differences between the two salts. One is the strength 
of the molecular dimerization, namely the dimerization 
in the I3 salt is larger than that in the SbFe salt. An¬ 
other difference is the transfers in the inter-stacking di¬ 
rection namely, the magnitudes as well as the sign of 
the inter-stacking transfers are different between the two 
salts, that is in c(g')-direction in the I3 (SbFe) salt. 

To clarify the origin of the differences between the two 
salts, we consider the alignments of the donor molecules 
in the conducting c-a plane for the two salts. The con¬ 
ducting c-a plane for each salt is shown in Fig. HI We 
find that the tilting angle of the donor molecules from the 
c-axis is different between the two salts. In the I3 salts 
shown in Fig. | 4 ] (a), the tilting angle is larger than that 
in the SbFe salts shown in Fig. | 4 ] (b). The difference in 
the tilting angle gives rise to differences in both the mag¬ 
nitude and the sign of the main inter-stacking transfers, 
which is tc in the I3 salt shown in the lower panel of Fig. 
m (a), while they are tqi and tq2 in the SbFe salt shown 
in the lower panel of Fig. | 4 ] (b). 

Now, let us try to decompose these differences. We 
consider a case where we hypothetically eliminate the 
dimerization effect. Namely, we simplify the model by 
considering only the nearest neighbor transfer energies. 
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TABLE I: List of the transfer energies in the unit of eV for 
^-(BDA-TTP)2X. 


A 

I 3 

SbFe 

tpi (eV) 

-0.174 

-0.153 

tp2 

-0.102 

-0.126 

tql 

0.018 

-0.071 

tq2 

0.041 

-0.055 

tc 

0.062 

0.007 

t2c 

0.002 

0.005 

t2p 

0.006 

0.021 

ta 

-0.001 

0.003 

t2q 

0.004 

0.005 

tq3 

-0.012 

0.003 

tq4 

0.013 

0.006 

trl 

0.002 

0.014 

tr2 

0.009 

0.008 

tp2/tpl 

0.586 

0.824 

tql -rtq2 

-0.214 

0.452 

tpl tp2 

2tc 

-0.449 

-0.050 



FIG. 4: (color online) The lattice structures of (a) the I 3 and 
(b) the SbFe salt. Blue solid lines represent the tilting angles 
of the donor molecules, which is measured from the c-direction 
taken in the horizontal direction. In the lower panel, the el¬ 
lipses represent the donor molecules, the black solid lines show 
the intra-stacking transfer, and the red solid (black dotted) 
lines represent the main (not main) inter-stacking transfer. 


and replace the hopping in the p- and g^-directions by tak¬ 
ing their averages. The band structure of the simplified 
model is given by 

s (k) = 2tcCos{kc)2tp cos {ka)21 q cos {kc 

where the transfer energies are tp = —0.138 eV, tq = 
0.030 eV, tc = 0.062 eV for the I 3 salt, and tp = —0.140 
eV, tq = —0.063 eV, tc = 0.007 eV for the SbFe salt. 
Eliminating the dimerization effect enables us to take the 
unit cell reduced along the a-direction. By comparing the 
band structure of the simplified model, we can extract the 
difference caused by the inter-stacking transfer. 

We compare the band structure of the two salts in the 
{kx, ky) plane, where ky is taken in the molecular stack¬ 
ing direction and kx is taken in the direction of the main 



FIG. 5: (color online) The simplified model for (a) the I 3 
salt and (b) the SbFe salt, where the line width schematically 
represents the magnitude of the transfer energies, and the 
solid (dashed) line represents the negative (positive) value of 
the transfer energies, (c) The band structure of the simplified 
model by eliminating the dimerization for each salts. 


inter-stacking transfer, namely, c-direction in the I 3 salt 
(Fig. El (a)) and a /2 + c-direction in the SbFe salt as seen 
in Fig. [5](b). Also, we shift the wave-number by (tt, 0) for 
the SbFe salt considering the sign difference in the main 
transfer energies along the inter-stacking direction. By 
such a transformation, we find that the band structures 
between the I 3 salt and the modified SbFe salt become 
very similar as shown in Fig. [5] (c). Since the simplified 
model eliminates the dimerization effect, the difference in 
the original band structure between the two salts comes 
mainly from the dimerization^ and the differences coming 
from the inter-stacking transfer are not essential. 


C. Effect of electron correlation cooperating with 
dimerization 

A quarter-filled system effectively becomes a half-filled 
system by increasing the dimerization^, so that the elec¬ 
tron correlation is strengthened. Since we now know that 
the strength of the dimerization is the essential difference 
between the I 3 and SbFe salts, we expect that the differ¬ 
ence of the ground state physical properties between the 
two salts is caused by the strength of the electron corre¬ 
lation originating from the difference in the strength of 
the dimerization. 

The strength of the electron correlation can be mea¬ 
sured by calculating the spin susceptibility. We apply the 
TPSC scheme to the Hubbard model of the I 3 salt. From 
the first principles calculation of the I 3 salt, the band 
width W is about 0.77eV, so we take the on-site interac¬ 
tion 1/0 = 0.8eV as same as the band width. The bare on¬ 
site interaction Uq is estimated in the other strongly cor¬ 
related organic conductors applying the first-principles 
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calculatioi>22i^. Referring to them, the on-site interac¬ 
tion we taken is appropriate. 

Figure [ 6 ] (a) shows the temperature dependence of the 
local vertex of the spin part Usp and the critical on¬ 
site interaction of the magnetic order Usbw fhe left 
scale. Above the temperature T O.OOdeV, Usp is almost 
unchanged and t/sDW gradually decreases with lowering 
the temperature. Below T « O.OOdeV, U^p takes almost 
the same value, but somewhat smaller value than t/sDW, 
which can be understood that the magnetic ordering is 
developed with lowering the temperature. 

In the right scale of FigEl^a), we present the ratio 
^SDw/tAp as a function of T. The TPSC approach sat¬ 
isfies the Mermin-Wagner theorem so the true magnetic 
ordering does not occur in the present model, but we can 
regard the temperature at which the line extrapolating 
I/sDw/^sp from high temperature reaches unity as the 
magnetic critical temperature in the actual three dimen¬ 
sional system. We estimate the magnetic critical tem¬ 
perature to be about 0.0038 eV. Reflecting the tendency 
toward the magnetic ordering, Uch. quickly increases be¬ 
low T = 0.0038 eV as shown in Fig. [ 6 ] (b). We show 
the double occupancy (D) = {n^ni) as a function of T 
in Fig. [ 6 ] (c). Similarly to the local vertices, the double 
occupancy (D) also changes below T = 0.0038eV. De¬ 
creasing the temperature reduces the double occupancy, 
which means the tendency of the magnetic localization at 
each site. Figure [ 6 ] (d) shows the inverse of the maximum 
value of the spin susceptibility against T. As expected 
from Fig. [ 6 ] (a), the inverse of the spin susceptibility ex¬ 
trapolates to zero around T = 0.0038 eV. In fact, a very 
recent experiment observes a magnetic transition in the 
Mott insulating state of the Is salt at low temperature^. 
TPSC is not capable of directly describing the magnetic 
ordering of a Mott insulator, but the very fact that the 
material is a Mott insulator is consistent with our view 
that the electron correlation effect is strong due to the 
strong dimerization. 

Figures 0 (a) and (b) show the absolute value of the 
Green’s function \G\ and the spin susceptibility Xsp of the 
I 3 salt with Uo = 0.8eV and T = O.OOdeV. The absolute 
value of the Green’s function takes large values near the 
Fermi surface shown in Fig. [7| (a). The wave number at 
which the spin susceptibility is maximized corresponds 
to the nesting vector of the Fermi surface as seen in Fig. 
[3(b). As shown in Fig. [7](b), the maximum value of the 
spin susceptibility takes a large value since its tempera¬ 
ture is close to the critical temperature. 

To clarify the relation between the electron correlation 
and the dimerization, we measure the strength of the 
dimerization by the quantity tp 2 /tpi • When tp 2 /goes 
to unity, the dimerization decreases. If the decrease of 
the dimerization results in weakening the electron corre¬ 
lation, we expect (i) U^p gradually deviates from I/sdw, 
(ii) the double occupancy {D) becomes large, and (hi) 
the maximum value of the spin susceptibility decreases 
within the TPSC scheme. Furthermore, if the stronger 
electron correlation of the I 3 salt originates from the 



FIG. 6 : (color online) Temperature dependence of (a) Usp 
and PsDW for left scale and Fsdq/Fsp for right scale, (b) 
Fch, (c) double occupancy (D), and (d) maximum value of 
the diagonalized spin susceptibility for the model of the I 3 
salt. Blue dashed lines represents the line extrapolating the 
each values and black dotted lines are about T = 0.0038 eV. 



FIG. 7: (color online) (a) The absolute value of the Green’s 
function and (b) the diagonalized spin susceptibility for the 
model of the I3 salt at T = 0.004eV. 


stronger molecular dimerization, all the quantities should 
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approach the values close to those of the SbF q salt when 
the dimerization is reduced hypothetically in the model 
of the I 3 salt. 

Let us now investigate the relation between the elec¬ 
tron correlation and the dimerization. Figure [ 8 ] (a) shows 
the local vertex of the spin part Usp and the critical on¬ 
site interaction for the magnetic order Usbw as a function 
of tp 2 ltpi in the model of the I 3 salt, also shows them for 
the SbFe salt at the point corresponding tp^j^px^ where 
we take T = 0.004 eV. Decreasing the dimerization (in¬ 
creasing tp 2 /tpi)^ Usp gradually differs from t/sDW, which 
expects that increasing the tp 2 ltp\ suppresses the maxi¬ 
mum value of the spin susceptibility. In contrast to the 
temperature dependence, decreasing the dimerization in¬ 
creases Uch as seen in Fig. [ 8 ] (b), although Usp differs 
from the t/sDW- la Fig- El (c), the double occupancy 
{D) monotonically increases with decreasing the dimer¬ 
ization, which can be understood as the suppression of 
the magnetic localization. This tendency is confirmed by 
the deviation of Usp from I/sdw- Figure [ 8 ] (d) shows the 
maximum value of the spin susceptibility as a function 
of tp 2 ltpi. Decreasing the dimerization from the actual 
value of the I 3 salt quickly suppresses the spin suscepti¬ 
bility, and that of the I 3 salt takes almost the same value 
as that of the SbFe salt around the same strength of the 
dimerization. 

From the tp 2 ltp\ dependence in Fig. [ 8 l we can say 
that the electron correlation in the I 3 salt is stronger 
than in the SbFe salt due to the strong dimerization. 
We therefore conclude that the difference of the ground 
state between the two salts, namely, insulating for the 
I 3 salt and superconducting for the SbFe salt, originates 
from the strength of the dimerization, which affects the 
electron correlation. Applying the pressure to the I 3 salt 
reduces the dimerization, resulting in the metallicity, and 
hence the superconductivity appears. 


IV. CONCLUSION 

In the present study, we have performed first prin¬ 
ciples band calculations and have derived the effective 
tight-binding models of /3-(BDA-TTP)2l3 and /3-(BDA- 
TTP) 2 SbF 6 . The band structures and the Fermi surface 
between the I 3 and SbFe salts are apparently different al¬ 
though only the anion differs. The derived tight-binding 
models, which accurately reproduce the first principles 
band structures of the two salts, show that the differences 
between the two salts comes mainly from the strength of 
the dimerization. 

As for the effect of the electron correlation, we have 
presented the TPSC results for quantities such as the spin 
susceptibility in the Hubbard model for the two salts. 
The TPSC results show that the electron correlation be¬ 
comes stronger upon lowering the temperature and/or 
increasing the dimerization strength. Then, we have hy¬ 
pothetically reduced the strength of the dimerization in 
the I 3 salt to that of the SbFe salt, where all the calcu- 



tplItpX 

FIG. 8 : (color online) The strength of the dimerization tp^jtpx 
dependence of (a) Usp and Fsdw, (b) Uch, (c) double occu¬ 
pancy (D), and (d) maximum value of the diagonalized spin 
susceptibility for the model of the I 3 salt with Uo = 0.8eV and 
T = 0.004eV. The each values on the SbFe salt with same U 
and T are also shown in the corresponding figures. 


lated quantities tend to become similar to those of the 
SbFe salt. Thus, we conclude that the electron correla¬ 
tion in the I 3 salt is stronger than the SbFe salt due to 
the strong dimerization. The expected stronger correla¬ 
tion in the I 3 salt is at least qualitatively consistent with 
a recent experimental observation that the material is a 
Mott insulator, which is a hallmark of strong correlation, 
and exhibits a magnetic transition at low temperature^. 
Applying the pressure to the I 3 salt reduces the dimeriza¬ 
tion, which weakens the electron correlation, and hence 
the superconductivity appears as in the SbFe salt. 

In the present study, we have considered only the on¬ 
site (intra-molecular) electron-electron interaction. It re¬ 
mains an interesting future problem to study the effect 
of the off-site interactions. In fact, it has been known 
that in organic conductors having quarter-filled bands, 
the Mott insulating state often competes with the charge 
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ordering and/or charge-density-wave states^. It is an in¬ 
teresting issue to investigate how such interactions would 
affect the insulating properties as well as the mechanism 
of the superconductivity. 
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